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^ ■ Abstract 

o ■ 

^ \ We study the spontaneous CP breaking at finite temperature in the Higgs sector in the 

r-| ■ Minimal Supersymmetric Standard Model with a gauge singlet. We consider the contri- 

p i' 

I ' bution of the standard model particles and that of stops, charginos, neutralinos, charged 

D . and neutral Higgs boson to the one-loop effective potential. Plasma effects for all bosons 

^ ■ 

are also included. Assuming CP conservation at zero temperature, so that experimental 

•T— I , 

^ . constraints coming from, e.g., the electric dipole moment of the neutron are avoided, and 

d \ the electroweak phase transition to be of the first order and proceeding via bubble nucle- 

ation, we show that spontaneous CP breaking cannot occur inside the bubble mainly due 
to large effects coming from the Higgs sector. However, spontaneous CP breaking can be 
present in the region of interest for the generation of the baryon asymmetry, namely inside 
the bubble wall. The important presence of very tiny explicit CP violating phases is also 
commented. 
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1. Introduction 



It is well known that for generating the observed baryon asymmetry of the Universe 
(BAU) three basic ingredients are necessary [||]: baryonic violating interactions, departure 
from thermal equilibrium and C and CP violation. Even if these conditions can be esaily 
fulfilled in grand unified theories (GUT's) 0], severe complications are present. First of 
all, any non zero fermion number {B + L) created at some superheavy scale is almost 
completely arised by the anomalous electroweak {B + L) violating processes , which are 
in equilibrium down to a temperature of about 10^ GeV. Moreover, if {B — L) violating 
processes are present and in equilibrium at high temperatures, then also an eventual (B — L) 
component of the BAU vanishes before the onset of the electroweak era, unless severe 
constraints on the {B — L) violating couplings are imposed [Q or ad hoc ways out are 
found §. 

All these considerations make the possibility of generating the BAU during the elec- 
troweak phase transition (EWPT) very appealing 0. Here the difficulty is threefold. First 
of all, it is still an open question if the necessary amount of CP violation to produce enough 
baryon asymmetry is present in the standard model (SM) 0; secondly, it is not obvious 
that the phase transition is enough strongly first order to make the baryon number produc- 
tion effective; third, one must avoid the wiping out of the baryon asymmetry by anomalous 
processes, which leads to an upper bound on the mass of the SM Higgs boson already 
ruled out by LEP results 0. 

The situation does not appear much more appealing when the minimal supersymmetric 
extension of the standard model (MSSM) ||10| is considered. The MSSM contains two extra 
CP violating phases with respect to the SM. The requirement that these phases provide 
the necessary amount of CP violation for the generation of the BAU gives rise to additional 



strong constraints on the parameter space of the model |TT[. Indeed, the electric dipole 
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moment of the neutron must be larger than 10"^'' e-cm, while an improvement of the 
current experimental bound on it by one order of magnitude would constrain the lightest 
chargino and the lightest neutralino to be lighter than 88 and 44 GeV, respectively. As far 
as the spontaneous CP breaking (SCPB) in the MSSM is concerned, it can be triggered by 



radiative corrections |]12[ at zero temperature. As established on general grounds by Georgi 
and Pais [1^, it requires the existence of a pseudoscalar Higgs boson with a mass of a few 
GeV [|I4|, which has been ruled out by LEP [0. Nevertheless, it has been recently realized 
[0 that temperature effects can trigger SCPB in the MSSM at the critical temperature 
of the EWPT and that the right amount of the baryon asymmetry can be generated with 
the help of tiny CP violating phases giving rise to a neutron electron dipole moment well 
below its present experimental limit |]T^ . Even if this mechanism can work in a wide region 
of the parameter space without any fine tuning, it requires, as a general tendency, small 
values of the mass of the pseudoscalar Higgs boson and large values of tan/? = vijvx^ v\ 
and V2 being the vacuum expectation values (VEV's) of the two Higgs doublets present 
in the model, whereas recent results on the phase transition in the MSSM seem to 
point out towards the opposite direction in order not to wash out the generated baryon 
asymmetry. 

In the framework of supersymmetric grand unified theories, the MSSM is not the most 
general low energy manifestation of supersymmetric GUT's. Indeed, it is possible that at 
low energy the theory contains and additional gauge singlet field, the so-called next-to- 
minimal supersymmetric standard model (NMSSM) |19|, as predicted in many superstring 
models based on and ^[/(S) ® |^ GUT groups. Incidentally, the presence 

of an additional gauge singlet allows to get rid of the so called /i problem in the MSSM, 
namely the presence of an unknown supersymmetric mass term in its superpotential. 

The study of the EWPT in the NMSSM has been performed in ref. where it 

has been shown that the order of the transition is determined by the trilinear soft super- 
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symmetric breaking terms and that it is possible to preserve the baryon asymmetry from 
the wiping out of sphaleron interactions for masses of the hghtest scalar well above the 
experimental limit. 

In this papar we want to address the question of the SCPB in the NMSSM at finite 
temperature, which is one of the key ingredients to generate the BAU at the electroweak 
phase transition. Indeed, in Supergravity inspired models with canonical Kahler potentials 
it is expected that CP violation in the parameters of the Higgs potential is not very 
significant ||23| and the presence of only explicit CP violation is not enough to produce 
the BAU, so that any CP violation in the Higgs sector must be spontaneous. Moreover, 



it has been shown by Romao that the potential of the NMSSM at the tree level and 
at zero temperature has no CP violating minimum. The point is that at the CP violating 
extremum there is always a mode with negative mass squared. When one-loop corrections 
to the tree level potential at zero temperature are added, the CP violating saddle point 



is turned into an absolute minimum pSf, where some of the Higgs fields are necessarily 
very light since they correspond to the modes having negative mass squared when one- 
loop corrections are not present. Quantitatively, the presence of a CP violating minimum 
requires two neutral and one charged Higgs bosons to have masses smaller than ~ 100 



GeV, which rules out much of the parameter space p5[. The aim of this paper is to 
investigate whether finite temperature effects can trigger SCPB in the NMSSM in regions 
of the parameter space which at zero temperature correspond to non CP violating minima, 
i.e. to mass eigenvalues for the lightest Higgses not dangerously light. 

The paper is organized as follow. In Section 2 we discuss the SCPB in the Higgs sector 
of the NMSSM on very general grounds, i.e. using the most general gauge invariant po- 
tential containing two Higgs doublets and the gauge singlet. In the same Section we also 
show how to renormalize the parameters present in the model giving all the details in the 
Appendix. This general discussion will make easier to understand why at zero temperature 
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CP cannot be broken spontaneously at the tree level and why, including zero temperature 
one-loop corrections, SCPB can occur only in a very restricted region of the parameter 
space. In Section 3 we introduce the corrections coming from the finite temperature effec- 
tive potential. In Section 4 we present and discuss our results making use of what learned 
in Section 2. Finally, Section 5 contains our conclusions. 2. Spontaneous CP 
Violation in the NMSSM 

2.1. General Analysis 

The superpotential involving the superfields Hi, H2 and in the NMSSM is 

W = XH1H2N - hN"" + htQH^U", (1) 

o 

where the A^^ term is present to avoid a global U{1) symmetry corresponding to N ^ Ne^^ 
and H1H2 — HiH2G~^^, and Q and U'^ denote respectively the left-handed quark doublet 
and the (anti) right handed quark singlet of the third generation. Note that a Z3 symmetry 
under which any Higgs superfield transforms as — acj) with a'^ = 1 is still present. This 
is also the typical structure emerging from superstring inspired scenarios. This symmetry, 
when spontaneously broken in the vacuum, could create a serious cosmological problem 
due to the appearance of domain walls. However, it has been shown in ref. that 
nonrenormalizable terms like ~ N'^/M, M is some superheavy scale as the GUT or the 
Planck scale, would prevent the density of domain walls from becoming large enough to 
create any cosmological danger, while being too small to have any impact in the low energy 
phenomenology as well as in the following discussions. 
The tree level potential is given by 

Vtree = Vp + Vd + Vsoft, 
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Vf = lAI^ 



lA^l' {\H,f + \H2f) + \H1H2f] + \kf \N\ 



\k*HiH2N*^ + h.c. 
Vd = U9' + 9")(m'-\H2fY + lg' HIH2 



2 



where 



2^ 

Ko/i = ml\H,f + ml\H2\^ + m%\Nf 

- {XAxHiH2N + h.c.)-(^hAkN^ + h.c.y (2) 



H.^l^l), (3) 



and g and g' are the gauge couphngs of SU{2)i and f/(l)y, respectively. 

Redifining the global phases of H2 and A^, it can be shown that all the parameters 
in eq. (2) can be made real, except the ratio r = A\/Ak. We assume this parameter to 
be real ||2^. Indeed, in Supergravity models with a canonical Kahler potential it turns 
out that r = 1 and then real at the unification scale, r can develop a phase through the 
renormalization group equations thanks to complex gaugino masses, but this effect is small 
due to the constraints on the gaugino phase coming from the electric dipole momentQ. 

If we define 

{Hi) = vi e''\ {Hi) = V2 e''\ {N) = x (4) 

and 

?.e^ = 2ip^, e, + e2 + e, = 2^^, e, + 62 - 29^ = 2ip, - 2^,, (5) 

we can write the most general gauge invariant (under SU{2)l ® U{1)y ® Z^) potential in 
the vacuum 

{V) = Divt + D2vt + D^vlvl + D^vlx" + D^vlx" + D^x^ 

+ DiViV2X^ COS {2ipi — 2(p3) + D^m\v\ + Dgm\v2 + -Diom^x^ 

+ Diif it;2X cos (2v9i) + Di2x'^ cos (2(/}3) , (6) 



^Note however that the presence of even very smaU phases for r is crucial for the generation of the 
BAU, see fllTll and befow. 



where it is easy to derive the D coefficients in the case of the tree level potential just 
comparing eqs. (2) and (6). 

The potential which depends on the phases 



V (v?!, v^s) = D'jViV2X^ cos {2(^1 - 2v?3) + -Diif 1^2^; cos {2(pi) + Di2X^ cos (2793) , (7) 



is minimized by the following values of ipi and <^3 ||25, ^ 

cos (2y9i) 
COS (2v?3) 
COS {2ipi — 2ip^) 
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where 



Ai2 = D11V1V2X, , Ai3 = DyViV2X^ , A23 = Di2X^, 

if |74i2v4i3|, 1^12^23! and |y4i3y423| form a triangle and if 

^12^13 



A 



> 0. 
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Defining now the new variable 



S = + t>2, A = vl — V. 



2! 



(9) 



(10) 



'IV. 



one can show that the extremum in the full parameter space (E, A, cp^) is a minimum 
if the following 3x3 matrix is definite positive 



where 



2i?4 Bq By 
B=\ Be 2B, Bs 
Bj Bs 2Bq 



D1 + D2 1^ , IDjDii 

IDA = U-X H — . 

^ 4 4 8 L)i2 



(12) 



6 





4- Do 
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Di 


-D2 




2 


D4 






2 




+ D5 



^ ^ 4 ^ 8 D12 



^6 = 

= 



2 ' 



B, = i^g-l:^. (13) 

For the tree level case, we have the particular values for the B parameters 

= A;2fi + iy (14) 



3r 

It is now easy to understand why CP cannot be broken spontaneously at the tree level 



Indeed, one has to satisfy the following system of conditions 



> =^ r < plus < 



123 



> ^ r < -| or r > 0, 
55 > ^ r > -|, (15) 
455^9 - 5| > ^ (1 + 3r)V3r > AVA;^ > 0, 



which does not have any solutions in the r space. The key point here is the following: once 
the condition given in eq. (10) is satisfied, which assures that the nonvanishing phases 
in eq. (8) correspond to a minimum, then the B matrix is never definite positive at the 
tree level. When one-loop corrections are considered, while the parameters involving the 
condition (10) receive small corrections because the latter are always proportional to A 
or k which are taken to be small (see below), among the elements of B, B^ may receive 
very large corrections from the top-stop sector and the Higgs sector (see the discussion of 
Section 4). It is then clear that, to have SCPB, the correction to the tree level potential 



should act predominantly on the B5 term to drive it and {AB^Bg — to positive values. 
This fact will be crucial for the following and expecially for the discussion about the finite 
temperature corrections to the tree level potential. 

2.2. How to Renormalize the Tree Level Potential 

In this Subsection we want to give some details on how to calculate the corrections to 
the tree level coefficients of the most general potential, eq. (6). A complete full analysis 
will be given in the Appendix. Note also that the following discussion can be extended to 
any model and to any potential. 

The one-loop effective potential can be written as 

Veff^Vtree+Vi, (16) 

where Vi H2, N) can indicate either the one-loop correction to the tree level potential 
at zero temperature or the one-loop correction at finite temperature. The SU (2) l'S'U (1)^0 
Z3 invariant and independent operators in the NMSSM are the dimension two operators 
\Hif, \H2f and \Nf and the dimension three operators 2 Re{N^), 2 Re{HiH2N*^) and 
2 Re{HiH2N). Let's collectively indicate them with £j and e^, respectively. The effective 
potential Vg// can be then written as follows 

i i<j i 

where Oj, Ojj and hi refer to dimension two, four and three operators, respectively. 

Since Vi receives its contributions from all the particles whose masses are field dependent 
and therefore is a function of such masses, one can easily calculate the coefficients of each 
gauge invariant operator in the following way (formulas for the dimension two operators 
are analogous to those for the dimension three operators) 
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dimension three operators: 2 Re{N^), 2 Re{HiH2N*'^) and 2 Re{HiH2N) 



dti 



dVt 



tree 



+ E 



dVi dml 



drni 9e,- 



dimension four operators: |-f/^i|^, |-f^2|^5 1^1^; l-^il^ l-^sl^, l-^il^ 1^1^ s-^d \H2f \N\' 



eff 



dsidej 

+ E 



tree 



dsidej 

1 2 2 



91^1 dml dml dVi d ml 



dml dsidsj 



In the above expressions the subscript indicate that all the derivatives must be evaluated 
for vanishing fields and the sum over the index a includes all the particles contributing 
to Vi. The complete details on how to calculate dml/dsi, dm^/dei, {d'^m\/ dsidsj) and 
{d'^ml/deidej) are given in the Appendix. 

Let's now consider the zero temperature contribution to Vi, V^^^. In the DR scheme 
of renormalization, it reads 



T=0 



1 



647r2 



Str <^ 



In 



(17) 



where Ai'^{(f)), with </> = {Hl,H2,N), is the field dependent squared mass matrix, the 
supertrace Str properly counts the degree of freedom, Q is the renormalization point and 
the dependence in eq. (17) is compensated by that of the renormalized parameters, so 
that the full effective potential is independent of up to the next-to-leading order. 

If we consider in V^^^ only the largest contributions [^, namely those coming from 
the top and the two stops with masses 



m. 



m 



thtR 



Ml + hi 



where we have taken a common soft mass term Ms for the left and right stop, one can 
show that the parameter receives, at the minimum, a large and positive correction 
proportional to h'l In (M|/mj) |2^. This is the reason why, when large zero temperature 



one-loop corrections to Vtree coming from the top-stop sector are considered, one can achieve 



SCPB only for large values of Ms This result does not contradict the Georgi-Pais 

theorem [|l^ which is based on the assumption that one-loop corrections to Vtree are small. 
Indeed, here one-loop corrections are large enough to turn the CP violating saddle point 
into a minimum. Note also that in the limit A ^ 0, the singlet tends to decouple from 
the {Hi, H2) sector and a CP violating minimum is found for {ipi, Lp^) = (±7r/4, ±7r/2) for 
Ak < and Ak > 0, respectively. In such Peccei-Quinn like symmetry is present 

in the superpotential involving Hi and H2 and it corresponds to the flat direction and 
then to the massless (at the tree level) pseudoscalar boson predicted by the Georgi-Pais 
theorem. This is the reason why in ref. |^ CP violating minima are found only for small 
values of A. 

Let's now come to the renormalization of the D coefficients at finite temperature. 
3. One-loop Corrections at Finite Temperature 



The one-loop effective potential at finite temperature is given by |27 



Ve 



Vtree + V^=' + V^^\ 



fer 



bos 



E 



n. 



48 647r2 
m,(0)2T2 m,(0)4 / 



In 



24 



647r2 



In 



AfT^ 2 
AbT^^ 2 



T 

12^ 



(mi(0)2 + n, 



3/2' 



where Ah = l6Af = 16 (3/2 — 27^;), •je being the Euler constant, and 



n. 



iJib) 



counts 



the effective fermionic (bosonic) degrees of freedom. Note that in the bosonic part Ilj 



denotes the thermal polarization mass for bosons contributing to the Debye mas [^. It 
arises when one resums at least the leading infrared-dominated higher-loop contributions 
to associated to the so called daisy diagrams pO) whose inclusions amounts to a 

resummation to all orders in a ~ {g'^ /27t){T^ /rn^), where g and m denote the generic 
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gauge coupling and mass respectively, and to neglect subleading contributions controlled 
by the parameter (3 ~ {g^ /27i)(T/m). 

Note also that in eq. (18) we have made use of the high temperature expansion for V^ff 



which turns to be a very good approximation for mif{p){(j)) < 1.6 (2.2) T ||30 |. 

Accordingly to what explained in the previous Subsection, the corrections to the coef- 



ficients and aij result as follow 



dimension three operators: 2 Re{N^), 2 Re{HiH2N*'^) and 2 Re{HiH2N) 



h, 



+ 



dei 
327r2 
327r2 



tree 



dei 



In 
In 



a,/er 

3^ 



48 



AfT^ 2^ 
A,T^ ^ 2 . 



. a,bos 

T 

~ 8^ 



dei 



24 



1/2 



(19) 



J 



dimension four operators: \Hi\ , \H2\ , |A^| , l-f^il 1-^21 , \N\ and \H2\ |A^|' 



eff 



d'V, 



tree 



+ 



deidsj 

a,fer 

d'^ml 
dsidsj 



dsidsj 



1 drn^ dml 
327r2 dsi dsj 



In 



A^T^ + 2, 



+ J2 ^a-i 



a,bos 



48 ^ 32^ 
[ (9m2 (9m^ 
1 dsi de 



1 



327r2 



In 



3 



T 



IGtt 



-1/2 



dsidSj 



+ 



T2 m2 



24 327r2 \ 



T 

8^ 



1/2 



(20) 



Working in the 't Hooft-Landau gauge and in the DR scheme, in the bosonic part we 
have summed over gauge bosons, stops, neutral and charged Higgs scalars, whereas in the 
fermionic part we have summed over top, neutralinos and charginos. Again we refer the 
reader to the Appendix for the complete analysis of all the technical points. 

The coefficients hi and aij given in the above eqs. (19) and (20) are the ones which 
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determine the D coefficients in the eq. (6) and then determine whether CP may be spon- 
taneously broken at finite temperature. 

In the next Section we shall focus on the results about the SCPB obtained when finite 
temperature corrections are taken into account showing that the picture can be consider- 
ably difi^erent from the case in which only T — Q corrections are considered. In particular, 
we shall address the question whether it is possible to have SCPB at finite temperaure in 
regions of the parameter space which correspond to CP conserving minima at T = 0, thus 
avoiding any limits coming from having very light Higgs bosons or electric dipole moment 
of the electron and neutron close to their experimental bounds. 

4. Results and Discussion 



As we have stressed in the Subsection 2.1, to have SCPB the corrections to the tree 
level potential should act predominantly on the term, see eq. (15). As a matter of 
fact, since < (when r < —1/3) at the tree level, the corrections to V^ff should be 
large enough to trigger positive values for S5. This is what happens in the case of T = 
top-stop large one-loop corrections, as pointed out in the Subsection 2.2. 

In the case of finite temperature corrections to Vtree-i we have proceeded as follows. First 
of all, we have fixed the set of the parameters given by (A, /c, Mg, m\, tan /3, At), where 
Mu is the soft breaking mass for the stop right and At is the trilinear soft breaking mass 
relative to the htQH2U'^ term in the superpotential. Imposing the minimization conditions 
at T = 

^ + ^ = 0, {<p, = v„<P, = v,,<P, = x), (21) 

allows one to express the parameters (^4^, Ax, m|) as functions of the already fixed param- 
eters and of two free parameters which we have chosen to be x and Mg, the latter being 
the soft breaking mass of the stop left. In the minimization conditions we have included 
only the relevant contributions to V^^'^, namely those coming from the top-stop sector. 
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The corresponding bottom contributions are negligible for tan f3 ^20 and the gauge sector 
does not play any important role. This also holds for the extended Higgs sector since the 
Yukawa couplings A and k are taken to be small as suggested by the requirement that they 
remain in the perturbative regime up to a large scale (say 10^^ GeV) ||19| or to forbid the 



breaking of electromagnetism |^T|. Note, however, that even if the Higgs sector does not 
play any role in the T = corrections, it will be fundamental in the finite temperature 
one-loop contributions as we shall explain later. 

In Figs, la) and lb) are displayed the plots of the experimentally allowed regions in 
the [Mq, x) plane for fixed values of the other parameters^. 

The first constraint is that the lightest Higgs CP even particle has not been produced 
in the decay Z° — > + h^. This gives the conservative bound rrifiO > 60 GeV dashed 
line (the exact bound depends on the coupling of to Z^, -R^o^o/jO, and then it is weaker, 
m/jO ^-R|o20/iO 60 GeV). In addition, the lightest pseudoscalar should be heavier than ~ 
20 GeV 1^ and the dashed-dot line corresponds to m^o = 40 GeV. Note that the allowed 
regions correspond to T = CP conserving minima. 

When finite temperature corrections are added to Vtree, one can define the critical 
temperature Tc as the value of T at which the origin of the field space becomes a saddle 
point for the effective potential. 



Det 



Mr^\T^)]^^, = 0, (22) 



where the effective mass matrix is given by the second derivatives of the full one-loop finite 
temperature potential with respect to the scalar fields. In the origin of the field space and 



^Very recently Ellwanger, R. de Traunbenberg and Savoy have scanned the complete parameter 
space searching for the allowed region compatible with different constraints. In particular, they have 
found that only very large values of x are permitted, x > 800 GeV. On the other hand this result is based 
on different assumptions, e.g. universality in the soft supersymmetry breaking at the GUT scale. As a 
consequence, we believe that it is safe to consider also smaller values for x. 
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in the basis {ReHl/ ^,ReH^/ ^,ReN / it reads 



MY^\T,) = Diag (m?,m^,m^y, 



+ -{X' + e) T\ (23) 



The critical temperature is then given by the highest temperature for which one of the rnf 
{i = 1,2, N) vanishes. Obviously, the effective potential becomes flat at the origin only 
along directions corresponding to negative soft masses. Due to the heavy top ml can 
run to negative values at low energy, whereas, for small A and k, both ml and remain 
positive. As a consequence, the EWPT is expected to occur first along the H2 direction 
and immediately afterwards nonvanishing VEV's for and N are driven by the Ax term 



In Figs, la) and lb) we have fixed the temperature at = 150 GeV and shown the 
curve corresponding to ml{T^ = 0, solid line. The points in the {Mg,x) plane lying on 
the curve are then the points for which Tc = 150 GeV. Indeed, since the EWPT is known 
to be first order, it occurs when ml is still positive, i.e. at a temperature T* higher than 
Tc'. since all the points in the {Mg,x) plane below the solid line correspond to ml > 0, 
then they correspond to the region of the parameter space where the EWPT occurs at 
temperatures smaller or equal to 150 GeV. 

We want to point out that the contribution from the Higgs sector, where plasma effects 
are crucial to make the effective potential real at the origin, does have an infrared sin- 



gularity proportional to {T /m2), responsible for the failure of perturbative expansion |27 
for values of {Mq,x) such that ml = 0. The curve fnl = is then representing an upper 
bound contour, for fixed T, in the {Mg, x) plane which severly contrains the region of the 
parameters. We have also checked that in the experimentally allowed regions below the 
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curve ml = and for the choice of (Mg, x) adopted for figures 2a) and 2b) (see below), the 
theory remais perturbative in the sense that the perturbative expansion (3 ~ g'^{T/2TT rn2) 
remains smaller than 1. 

Once the temperature has fallen down to the tunneling temperature T*, the EWPT 
proceeds via bubble nucleation. In the centre of the bubbles, to avoid the wiping out 
of the baryon asymmetry by anomalous interactions, we have to impose that v{T*) = 
^Jvl(T*) + f|(T*) > T*, which is easily satisfied in the model under consideration p2| . 
Even if the SCPB in the bubble is not interesting as far as the generation of the baryon 
asymmetry is concerned^, it is however interesting per se ' asking whether the system can 
tunnel first to a CP violating minimum and then, as the temperature decreases, reach the 
CP conserving minimum which represents our initial condition at T = 0. 

4.1. SCPB inside the Bubble 

Scanning several sets of points in the allowed region in the [Mq, x) plane, we have 
numerically minimized the effective potential and looked for portions in the plane where 
conditions (15) could be satisfied. This research has turned out to be fruitless in the sense 
that very small and, in practice, negligible, allowed regions in the [Mg, x) plane have shown 
up. Even if this result might be quite surprising at a first sight, it can be explained as the 
effect of different phenomena happening at finite temperature. 

We know that one- loop corrections to Ve// should trigger positive values for in order 
to have SCPB. First of all, we have seen that at zero temperature the top-stop sector gives 
a large positive contribution to proportional to ln(M|/m^) for large Ms- At finite 
temperature, the corresponding contribution is reduced to hf \n{Ab/Af). 

Secondly, the Higgs sector plays a crucial role in determining the sign of B^. Us- 

In the usual scenarios d for the generation of BAU at the EWPT, the only source of baryon violations 
lies in sphaleron interactions whose rate is imposed to be small enough in the propagating bubble. 
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ing eqs. (6), (13), the formalism developed in the Appendix and taking into account 
that near the critical temperature m2 <^ rni,rn^, one can realizes that the corrections 
to i?5 are dominated by negative terms such as —{T/m2){dfn2/d\H^\^)'^, with i = 1,2 
and by —{T/m2)idfn2/d\Hff){dfn2/d\H2f), so that the largest Higgs contribution to 

i?5, AB^^max, is 

1 ^ W25 4 30 ,4 19 2 ,2\ r. 
Svr 7712 16 V 4 4 2 J 

Thus, -B5 receives a large and negative contribution from the one-loop corrections at finite 
temperature from the Higgs scalar sector. Since this contribution turns out to be, in 
absolute values, the largest among the different particles, this prevents to become 
positive and, consequently, the matrix B to be definite positive. This is the main reason 
why SCPB cannot occur in the centre of the propagating bubble. 



4.2. SCPB inside the Bubble Wall 



Fortunately, the situation changes when one moves away from the centre of the bubble 
towards the bubble wall. As a matter of fact, since the EWPT is of the first order, the 
temperature keeps constant until the Universe is in broken phase, whereas fi, V2 and x 
change their values from zero to fi(T*), V2{T*) and x{T*), respectively, when a bubble 
wall passes through a fixed point in space. Incidentally, we define the bubble wall as the 
region in which sphalerons are active, that is in which < v{T*) < T* and we remind the 
reader that, since here the BAU can be created through, e.g., the reflection baryogenesis 
mechanism [Q, it is just in this region that one needs SCPB. 

To describe what is happening inside the bubble wall one should solve a system of 
differential equations involving vi, V2, x, ipi and ips at T*. For instance, the equation for 
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(^3, should read 

d 
dz 



x\z,T*)-^^cp,{z,T*) 



4 dVeff , , 



where we are approximating the bubble wall to an infinite plane propagating perpendicu- 
larly to the z axis with a width L^. Nevertheless, one can make the analysis much simpler 
envisaging the following situation. In the unbroken phase and up to the edge of the ad- 
vancing bubble wall, at z = z^, f i and V2 are constantly equal to zero and then there is no 
SCPB; in the broken phase and up to the other edge of the bubble wall, at z — z^ + L^j, 
Vi and V2 are nearly constantly equal to their minimum values at the centre of the bubble 
and, as shown before, again no SCPB is occuring; finally from inside to outside of the 
bubble wall, the VEV's are decreasing from their values inside the bubble to zero. 

The key point now is that in the bubble wall, where equations like eq. (25) are valid, 
the VEV's Vi{z,T*), V2{z,T*) and x{z,T*) do not have to satisfy any longer the severe 
constraint Bq > coming from the requirment of a mass matrix with positive determinant, 
but one can simply impose that at the two edges of the bubble wall, where, e.g., eq. (25) 
reduces to dVeff / d(pi{z,T*) = 0, cos2ip3{zw + Lyj,T*) > 1 and cos 2933(2;^, T*) < -1. So 
doing, even if the solution of eq. (25) is not known, one is assured that maximal SCPB 
is occuring inside the bubble wall, that is an observer in a fixed point in space would 
experience a change Ay93 = 7r/2. Imposing cos2(^3(z^^, + L^,T*) > 1 and cos 2(^3(2;^, T*) < 
— 1 at the edges of the bubble wall is equivalent to impose 

0<t)_ <t;(T*) <{}(^^ + L^,T*), (26) 

where v± are found imposing cos2ip^{T*) = ±1, respectively, and are given by 



v± = 



%/2 



±^ x\z, T) + ^ x^z, n + ^ x\z, n 



1/4 



(27) 



I sin 2/^1^/^ 

At v{z,T*) = v± eq. (25) is then equivalent to dVeff /dips{z,T*) = and also eq. (10) 
must be satisfied. Naturally, similar relations are vahd for v'^{z,T*) and x{z,T*) once one 
imposes similar conditions on the angles 2(pi and 2{(pi — (p^). 
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In Figs. 2a) (2b)) we present the region in the {v{z,T*),x{z,T*)) plane corresponding 
to the conditions (26-27) for the choice x = 600 (800) GeV, niN = 10, 300 GeV and 
Mg = 200 (300) GeV. The points inside the regions I and II correspond to values of 
•v{z,T*) and x{z,T*) for which a maximal SCPB is occuring in the bubble wall, i.e, at 
least one of the angles 2^93, 2(pi and 2((/?i — ips) change by an amount vr inside the bubble 
wall. 



As pointed out in ref. [T^, the only presence of SCPB inside the bubble wall is not 
sufficient to create a net BAU. Indeed, if CP is spontaneously broken, any phase, let's call 
it 6, can take two opposite values, corresponding to two exactly degenerate vacua, since 
the effective potential depends only on cos 5. This degeneracy between the two vacua ±5 
would induce an equal number of nucleated bubbles carrying phases with opposite signs, 
which in turn generate baryon asymmetries of opposite signs and an overall BAU equal 
to zero when averaging over the entire volume of the Universe. However, the presence 
of a very small and explicit phase in the effective potential, for instance in r, of order of 
10"^ - 10-^ is sufficient to lift the de generacy leading to a difference between nucleation 
rates of the two kind of bubbles and to a baryon asymmetry of the right order of magnitude 



1^. This is due to the fact that, being the nucleation rate proportional to exp(— AF/T), 
where AF is the difference in free energy between the two vacua with phases ±5, it is 
very sensible to even small changes in AF. The smallness of the explicit phase does not 
change quantitatively our conclusions on SCPB and give rise to negligible contributions to 
the neutron electric dipole moment. 

5. Conclusions 



In the present paper we have investigated the possibility of SCPB at finite temperature 
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in the MSSM. After having performed a systematical analysis of the renormalization of the 
operators present in the potential, we have shown that SCPB at finite temperature cannot 
occur inside the propagating bubble walls which appear after the EWPT. This is due to 
the fact that the Higgs contribution turns out to be crucial and pushing towards the wrong 
direction. 

SCPB can occur inside the bubble wall, which is the interesting region as far as the 
generaton of the BAU is concerned, without any fine-tuning and in regions of the parameter 
space which corespond at T = to no SCPB, i.e. to regions where no constraints are 
present coming from experimental bounds on the electron and/or neutron electric dipole 
moment and from searching for the lightest Higgs bosons at LEP. 

Very small explicit phases are necessary to generate a nonvanishing BAU and in any case 
they give rise to an amount of CP violation which is well below the current experimental 
limits. 
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Appendix. 

In this Appendix we want to give a complete description of all the techniques necessary to 
calculate the corrections to the coefficients and a^j. Note that this description is quite 
general and does not depend at all on the particular model we are working with. 

As explained in the text, to calculate the renormalized and a^j coefficients, we need 
to calculate dml/dei and (d'^ml/deidej), ml being the field dependent mass eigenvalues of 
the particles contributing to the effective potential. Since are by definition the squared 
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mass eigenvalues of an A^xA^ hermitian mass matrix they obey the following equation 



N 



u 



= Det 



MU<I>) - ^a(0)J = E = 0. (A. 1) 

n=0 



Let's analyze first the case for a non degenerate mass matrix Taking the derivative 

of both sides of eq. (A.l) with respect to the gauge invariant and independent operator Sj, 
one can easily finds that (analogous expressions are valid for the dimension three operators 

(A. 2) 



2{n-l) ' 



while, taking the second derivatives with respect to the operators and Sj, one finds 



dSidSj 



N 



1-1) E 



v^AT 2(n-l) 
Z^n=0 TT'CnUla n=0 



n{n - l)c„ 



dml 9mg 2(n-2) 



+ m 



2n 



2(n-i) ( ^9ml , dcndml^ 



(A. 3) 



In the case of degenerate squared mass matrices, e.g. the 6x6 squared mass matrix for the 
neutral Higgs bosons (see below), one cannot use the equation 



dsi dml 9ei " dsi 







(A. 4) 



to invert in favour of [dm^/dsi) since now {dU/dm\) is vanishing. Nevertheless, for 
eigenvalues with degree of degeneracy equal to two, one can take the derivative of eq. (A.4) 
with respect to and S^m^/Sei^jthen obtains 



dml_ Etonm^("-^)fe 



(A. 5) 



and 



dsidsi 



2(n-2) 



N 

E 

n=0 L 



n(n — l)(n — 2)c, 



dml gmg ^2(„_3) 



+ nm 



2(n-l) 



- — h n(n — 



2(n-2) 



9c„, ^ dcn dml 



(A. 6) 
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What one really needs to calculate the corrections to the coefficients and a^-, which 
in their turn determine the possibility of having SCPB at finite temperature, are then 
the coefficients for each squared mass matrix whose eigenstates contribute to the effective 
potential. In the following we give all the details about the properties of the particles 
which must be included in Veff and from which the coefficients c„(0) can be extracted. 

• Top: Tit = -12, m2 = h^^ \H^\'^. 

• Stop: n^^ = ni_^ = 6; the squared mass matrix is 



Mt 



m 



LL 



m 



LR 



(A. 7) 



where 



m 



LL 



m 



RR 



M? + Us + h 



Hi 



2 

\12 4 , 



Hi 



Mf + + k 
ht (AtHl + XN*Hf 



9 



12 



Hi 



Hk 



(A. 



and 118 



2r7-i2 



9 



9tT 



(A. 9) 



are the thermal polarization squared masses for the squarks calculated in the limit of heavy 
gluinos. 

• Gauge bosons: nw± = 6, nw± = 2, n^o = 3, n^o = n^^ = 1, 



m 



9_ 
2 



9^ + 9 



Hi 

12 



+ 



H' 



Hi 



Hi 



m|o 



± 



1 



g2 _ g>2 



1 2 



Hi\' + \Hin+u^±-u 



r0|2 



Bl 



+ 



99_ 
2 



(A. 10) 



{\Hlf + \Hl\ 



n 2 
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where ITH 



5 



47 
10 



/rr2 



(A. li: 



Note that only the longitudinal components of the gauge bosons do have a Debye mass 
proportional to T^. 

• Charginos: = 2; in the basis {—iX~^, —iX~, ipfj-^) the mass matrix is 



f M2 gHf \ 

gHf \N 

Ms gHf 

V gHf XN J 



(A. 12) 



where M2 is the gaugino soft mass for the gauge group SU{2)l. 

• Neutralinos: n^o = 2; in the basis {W^, 5°, H^, H2, N) the mass matrix is 



/ 


M2 


- ^Hf 


a ij*o 

2-"2 





\ 







Ml 


a' u*o 


9' C7*0 
T-"2 







g zj*o 


9' ZJ*0 





AA^ 


XH^ 




a ij*o 


g' f7*o 
T-"2 


XN 





XH^ 


V 








XHI 




-2kN J 



(A. 13) 



where Mi is the gaugino soft mass for the gauge group U{1)y- 

• Charged Higgs scalars: ?2/^± = 2; in the basis {Hi,H2*) the squared mass matrix 
reads 



_ff±,ii 



ml + X'\Nf + Ug' + g'') (\h',\\\h'2 



H±,22 



4 

ml + X'\N\'-\{g- + g'^ 



Hi 



+ 



Hi 



V — Hi 
2 ^ 



A^^±,i2 = -M|^±,2i = AAaAT + XW + _ ^2^ EfHf. (A. 14) 

• Neutral Higgs scalars: rij^o = 1; in the basis {H^, Hf, H^, Hf, N, N*) the squared 
mass matrix is given by 



-^^0,11 = m^ + A'lA^l' + A 



Hi 



2 1 

'"4 



9' + 9" 



H'2 



Mlo^,2 = Mj^o,2i = Mlo^,, = l{g' + g")HlHl 



22 



•^H0,13 



•^^0,41 — -^^0,23 



(9' + 9") 



ttO tt*0 
-0 ttO \ 7. At2 



■M|;o,5i = ■M|;o,26 = ^'H^N* - XkH*"N, 



Mh022 


= m^ + A^lTVl^ + A^ 






Hi 




n 2 




= ml + X^\Nf + X' 


Hi 




Hi 




n 2 


•^H0,35 


= Mh0^53 = -^1^0 _46 


= A^i/^^* - XkHfN, 






•^H0,36 




= X^Hl - XAxHf, 










= m2 + A2|7V|2 + A2 


Hi 




Hi 




n 2 



+ A^ 



Hi 



+ A^ 



^2° 



+ 2A;2|iV|% 



M^o,56 = M%^Q^^2eN'' + 2XkHlHl-^kAkN\ 



•^^0,66 



mi + A' 



(-AT 



+ A^ 



Hi 



(A. 15) 
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Figure Caption 

Fig. la) For the following values of the parameters: A = 0.2, k = 0.4, ht = 1, = Mg + 10 
GeV, mjv = 10, 300 GeV, mi = 250 GeV, tan p = 2.5, At = 10 GeV, Mi = M2 = 150 
GeV and T — 150 GeV, and in the plane (x, M,), the solid hne corresponds to 
the curve m2{T) = 0, the dot-dashed line corresponds to m^o = 40 GeV and the 
dashed line to the conservative bound rriho^ 60 GeV. Regions I and II correspond to 
mjv = 10, 300 GeV, respectively. 

Fig lb) The same as in Fig. la) but for A = 0.3, k = 0.5, tan/5 = 1.2, mi = 150 GeV, 
Ml = M2 = 250 GeV. 

Fig. 2a) Values in the plane {v{z,T*), x{z,T*)) for which a maximal SCPB occurs. Region 
I and II correspond to the choice of the parameters given in Fig. la) with x — 600 
GeV and M, = 200 GeV. 

Fig. 2b) Values in the plane {v{z,T*),x{z,T*)) for which a maximal SCPB occurs. Region 
1 and II correspond to the choice of the parameters given in Fig. lb) with x = 800 
GeV and Mg = 300 GeV. 
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